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Abstract 

Propagating, solitary magnetic wave solutions of the Landau-Lifshitz equation with uniaxial, easy-axis anisotropy in thin (two- 
dimensional) magnetic films are investigated. These localized, nontopological wave structures, parametrized by their precessional 
frequency and propagation speed, extend the stationary, coherently precessing "magnon droplet" to the moving frame, a non-trivial 
generalization due to the lack of Galilean invariance. Propagating droplets move on a spin wave background with a nonlinear 
droplet dispersion relation that yields a limited range of allowable droplet speeds and frequencies. An iterative numerical technique 
is used to compute the propagating droplet's structure and properties. The results agree with previous asymptotic calculations in 
the weakly nonlinear regime. Furthermore, an analytical criterion for the droplet's orbital stability is confirmed. Time-dependent 
numerical simulations further verify the propagating droplet's robustness to perturbation when its frequency and speed lie within 
the allowable range. 
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1. Introduction 

Magnetic materials yield a rich variety of intrigu- 
ing nonlinear wave phenomena. Recent theoretical 
and experimental developments [1 4] have enabled 
the controlled manipulation of magnetic moments on 
the nanometer length scale, the magnetic exchange 
length, thereby generating further interest in the field 
of nanomagnetism [5]. Dynamic, strongly nonlinear, 
localized wave structures varying on the exchange 
length scale have been studied theoretically for some 
time [6] (see the exhaustive review [7]). It was re- 
cently proposed that one member of this family of 
solitary waves, the two-dimensional (2D) nontopo- 
logical droplet, could be experimentally realized in a 
spin torque nanocontact system that delivers a local- 
ized torque to a ferromagnetic thin film, balancing 
the inherent material damping [8] . The so-called dis- 
sipative droplet can, under certain conditions, expe- 
rience a drift instability, propagating for some time 
before eventually succumbing to magnetic damping. 
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A natural question arises, then, to inquire into the 

existence, stability, and properties of propagating, 
nontopological solitary waves in the underlying con- 
servative model, the Landau-Lifshitz (LL) equation. 

Previous studies of 2D propagating droplets in- 
clude numerical computations in the isotropic [9] and 
anisotropic [10] regimes as well as weakly nonlinear 
asymptotics [11]. In [10], a numerical study of topo- 
logical and nontopological localized solutions in 2D 
ferromagnets with easy-axis anisotropy was under- 
taken. Stationary droplets were constructed numer- 
ically which were then made to propagate in time- 
dependent simulations by the imposition of an ini- 
tial, unidirectional phase gradient. The following ap- 
proximate properties were observed: the excitation, 
accompanied by spin wave radiation, moved coher- 
ently with constant speed, a single frequency of pre- 
cession, and a spatially dependent phase lag. In this 
work, we demonstrate that the excitations observed 
in [10] and the moving, coherent structures in [8] can 
be identified as propagating, nontopological droplet 
solitary waves which we construct directly by study- 
ing a nonlinear eigenvalue problem. 
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We undertake a careful mimerical and asymptotic 
study of propagating 2D nontopological droplets, 
from here on in termed droplets. A two-parameter 
family of localized solitary wave solutions to the LL 
equation with easy-axis anisotropy parametrized by 
speed V and rest frequency w is constructed. In order 
to avoid a linear resonance, we arc led to postulate 
the droplet's existence and stability for speeds and 
frequencies in the restricted range 



y2 

l-w - — >0, 



similar to the existence conditions for moving ID 
droplets [7,12]. The droplet's structure is studied an- 
alytically in the large amplitude, small speed regime 
and in the weakly nonlinear regime following [11]. A 
nonlinear dispersion relation for a background spin 
wave which modulates the propagating droplet natu- 
rally arises. It degenerates to the linear dispersion re- 
lation for exchange spin waves in the small amplitude 
regime where the Nonlinear Schrodinger equation ap- 
proximately governs the dynamics. An iterative nu- 
merical method is employed to compute the droplet's 
structure in the comoving frame by solving an ap- 
propriate boundary value problem. These results are 
used to compute the amplitude, asymmetry, energy, 
momentum, and more for a wide range of propagat- 
ing droplets which agree with the asymptotic calcula- 
tions. Droplet orbital stability is confirmed by numer- 
ical verification of an analytical, Jacobian condition. 
Time-dependent numerical simulations with initially 
perturbed propagating droplets further demonstrate 
the robustness of droplet propagation. 

The outline of this work is as follows. In Sec. 2, 
several equivalent forms of the LL equation are intro- 
duced for later use in asymptotic and numerical anal- 
ysis. In addition, properties of the stationary droplet 
are reviewed. The propagating droplet problem for- 
mulation and some basic properties are discussed in 
Sec. 3. Section 4 presents the asymptotic results. In 
Sec. 5, numerical results addressing droplet proper- 
ties and stability are given. Further discussion and 
conclusions are related in Sec. 6 and 7. Finally, vali- 
dation of the numerical technique is presented in the 
Appendix. 



2. Model Equations and the Stationary 
Droplet 

The Landau-Lifshitz torque equation 



dM 



(la) 
(lb) 



OA / 2K 

;AM + ( i/o + -^M3 ]e3 + Hd, 



describes the dissipationless dynamics of the magne- 
tization M : X [-(5/2, S/2] x M+ M^S^ in an 
unbounded ferromagnetic film of thickness S exhibit- 
ing perpendicular, uniaxial anisotropy and a perpen- 
dicular applied magnetic field [13]. 

Relevant parameters for Eq. (1) are the gyromag- 
netic ratio (7), the free space permeability (/Uq), the 
exchange stiffness parameter (A), the perpendicular 
magnetic field amplitude {Hq), and the crystalline 
anisotropy constant (i^u)- 

The magnetization 



M = M,m ; 



(2) 



satisfies \M{x,t)\ = Mg, where Mg is the saturation 
magnetization, as can be verified by taking the dot 
product of (1) with M. The vectors ej, j = 1,2,3, 
are the standard Cartesian basis vectors for (see 
Fig. 1). The magnetostatic or dipolar field Hd is, 
in general, nonlocal, coupling the magnetization dy- 
namics to Maxwell's equations which, in the magne- 
tostatic regime, take the form 

V -ffd = -V • M, 

where M is extended to zero outside of M? x 
[—S/2, S/2]. We are interested in the thin film regime 
where 5 is sufficiently small. There are a number 
of different thin film scalings that depend upon the 
length scales inherent to the problem at hand. In 
the static case, some thin film limits have been rig- 
orously justified via F-convergence; see, e.g., [14,15] 
and references therein. Thin film magnetodynamics 
have also been rigorously studied [14,16]. In many 
thin film scalings, there are two key simplifications 
that can be made: the magnetization is approxi- 
mately uniform through the thickness of the film so 
that the problem is approximately 2D 



(3) 



and the dipolar field can be approximated by the 
local term 

Hd ^ -Maeg. (4) 

Heuristieally, the two-dimensionality of the problem 
is justified when the film thickness is much smaller 
than the transverse magnetic excitations of interest 
because the exchange energy for magnetization varia- 
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tions in the .T3 direction would be proliibitivc^ly large. 
In our case, this amounts to the requirement 

where L^x is the magnetic exchange length 



(5) 



1/2 



Q = 



typically on the order of several nanometers, and Q 
is the dimensionless quality factor 

assumed here to be greater than unity. Assuming that 
dx^M = 0, we can compute the 2D Fourier transform 
of the X3 averaged dipolar field (see e.g. [17]) 

k[k-M_L{k)]. 



{Hd + M.es) 



A;2 



■[i-r(M)] 



(6) 



where /(fc) = /^^ f{x±)e 
M_L = (Mi,M2), and 

f(M) = 



-(M3-M,)(fc)r(A;<5)e3, 

dXji_, X± = {xi,X2), 



- e~ 
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(7) 



The Fourier transforms are computed for rapidly de- 
caying functions. In this work, we study solutions 
that rapidly decay to M MgE^ with the decay 
length Lex/iQ — 1)^^^ (see non-dimensionalization 
below). Then the Fourier transforms in (6) are con- 
centrated on wavevectors k of order {Q — l)^^"^ /Lex 
or less. Under the assumption (5), we can approxi- 
mate (6), (7) as 



r(M) ~ 1 - -M = 1 + O U 



Hd = -Mgeg + 0{5 



hence (4) is justified. 

It is convenient to nondimensionalize Eq. (1) ac- 
cording to (2) and 



iO- 11^/2 

t' = |7|MoM,(Q-l)i, x'= ^^ ' X, 



ho 



Ho 



(8) 



Ms(Q-l)' 

so that, after dropping primes, the problem of study- 
ing localized solutions to Eq. (1) subject to (3) and 
(4) becomes 

dm 



at 

m : 



-m X [Am -I- (ma -I- ho) 63] 
X R+ ^ S^, lim m{x, t) = 63, t> 0. 

|a!|— >cx) 

(9) 



The nondimensionalization (8) assumes an easy-axis 
anisotropy so that the quality factor Q, measuring 
the relative strength of the anisotropy to the magne- 
tostatic energy, is larger than unity. When Q < 1, the 
anisotropy is of the easy-plane variety. Anisotropic 
ferromagnetic nanostructures in this thin film geom- 
etry can be fabricated, e.g., using Cobalt/Nickel mul- 
tilayers [5]. 

Transforming to the rotating frame 

m! = [mi cos {ho t) + m2 sin {ho t), 

— mi sin {ho t) + m2 cos {ho t), ms], 

takes Eq. (9) to the parameterless form (after drop- 
ping the prime) 

dm 



dt 



= -m X heflf, /ieflF = Am-l-mses. (10) 



Therefore, without loss of generality, we will consider 
localized solutions of Eq. (10). 

Equation (10) admits a family of localized, coher- 
ently precessing solutions [6,7]. The solutions that in- 
corporate a complete reversal of the magnetization, 
i.e. m(a;o, t) = and lim|a,|^oo m{x, t) — ±63 for 
one Xo e (magnetic vortices), are called topolog- 
ical and exhibit pinning of the soliton center, hence 
cannot propagate [18]. In contrast, solutions that do 
not exhibit magnetization reversal are termed non- 
topological and can propagate. Such solutions are the 
focus of this work. 

Sufficiently localized solutions of Eq. (10) conserve 
the magnetic energy 



5[m; 



|Vm| V (1 - m^) 



da; . 



(11) 



which incorporates contributions due to exchange 
|Vm,p and anisotropy (1 — m|). Integration is always 
assumed to be taken over the plane unless other- 
wise noted. The effective field in Eq. (10) is derived 
from the magnetic energy via 

h 

dm 

Other conserved quantities include the spin density 
J\f[m] = y^(l - ^3) dx, (12) 



and the momentum 



V[m] 



m2Vmi — miVm2 



dx. 



(13) 



1 -I-TO3 

The quantity Af can be interpreted as the mean num- 
ber of spin deviations from the uniform state 63 in a 
localized magnetic excitation. In the quantum limit 
it takes integer values and corresponds to the num- 
ber of magnons in an excited ferromagnet [7] . 
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Fig. 1. Schematic of tiie coordinate systems for a planar fer- 
romagnet consisting of the domain of m (R-^ in polar coordi- 
nates with radius p, azimuthal angle ip) and the range of m 
(unit sphere with polar angle 0, azimuthal angle 

Note that the momentum as given in (13) is am- 
biguous when the magnetization involves topologi- 
cal structures. The ambiguity has been addressed by 
considering moments of the topological vorticity [18] . 
Since we are considering nontopological structures, 
the momentum in (13) is equivalent to the more gen- 
eral formulation of [18] . 

Small amplitude exchange spin waves deviating 
from the background state 63 in the form 

mi = ecos(fe • X — uji^t), 

m2 = —esin{k ■ X — uji^t)], 7713 ^ 1, jej ^ 1 , 

with wavenumber k and frequency wl satisfy the lin- 
ear dispersion relation 

a;L(fe) = WL(fc) = l-|-fc^ k = \k\. (14) 

2.1. Alternate Forms of the Landau- Lifshitz 

Equation 

Two alternative forms for Eq. (10) will be used in 
this work. For the first, we use the spherical basis (see 
Fig. 1) 

m = [cos $ sin O, sin $ sin ©, cos 0] , 

where $ € [0, 2-n) (mod 2tt) and 9 G [0, tt] are the 
azimuthal and polar angles, respectively. Inserting 
this expression into Eq. (10) gives 



~dt 
where 



F[e,$] = 



F[e,$], 

V- (sin^eV^) 



sine— =G[e,$], (15a) 



sinG 



(15b) 



G[e, $] = - sin(2e) (] V^l^ + 1)-AQ. 



We will see that Eq. (15) is in a useful form to study 
stationary and slowly moving droplets as well as to 
formulate an energy minimization problem for prop- 



agating droplets. The energy, spin density, and mo- 
mentum in these variables take the form 



f[e,$] 



jVej^-hsin^e 



(l + |V$|^) 



M[Q\ = j (1-cose) dec, 
>[e,$]= j (cosO-l) V$da;. 



da; , 

(16a) 
(16b) 

(16c) 



Note that (15) is Hamilton's equations in the form 
__ S£ dcosQ _ 5£ 

Ik ~ ~(5cose' dt ~ (5$' 

where $ and cos 9 are the canonically conjugate vari- 
ables for the Hamiltonian system. 

Another representation of (10) can be had by per- 
forming the stereographic projection 

nil + ini2 e** sin 9 



w 



1 -|- ms 1 -|- cos 9 
Then (10) becomes [19] 

2w*V«; • Ww + w{l - \w\'^) 



iwt — Aw 



1 



w 



(17) 



(18) 



where the appended asterisk denotes complex conju- 
gation. This form of the LL equation is amenable to 
weakly nonlinear asymptotics (Sec. 4) and numerical 
analysis (Sec. 5). The conserved quantities are 

|2 , I ,2 



£[w] = 2 

Af[w] = 2 y - 



Im(w* Vw) 



da; , 



dx . 



(19a) 
(19b) 
(19c) 



In analogy with a Bose-Einstein condensate [20] , the 
momentum (19c) corresponds to a weighted phase 
gradient or a "magnetic superfluid" momentum. 

For the time-dependent initial value problem dis- 
cussed in Sec. 5.3, we will solve Eq. (10). Thus, we 
will make use of all three forms (10), (15), and (17) 
of the LL equation in our analysis. 

2.2. Stationary Droplet 

An azimutlially symmetric, stationary, localized 
droplet solution to the LL equation was studied in 
[6,7] . We note that the term magnon droplet was used 
because of its analogy with a localized "droplet" of a 
condensed, attractive one-dimensional Bose gas [7]. 
In this section, we briefly review the properties of the 
stationary droplet. 
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It is convenient to write the stationary droplet so- 
lution in a spherical basis satisfying Eq. (15). It takes 
the form 

Q{x,t) = Qo{p;uj) , ^{x,t)=wt, 
F[Qo,cot\=0, G[eo,ujt]=uj smQo, 



where p = ■\/xf~+~x| is the distance from the origin 
and the polar angle function Qq, parametrized by 
the droplet frequency w, is determined by solving the 
following nonlinear eigenvalue problem 



dp2 



1 d 
pdp 



dp 



(0;w) = 0, 



Qq + sin 60 cos 60 — w sin 60 = , 

(20a) 

lim eo{p;Lo) = 0. (20b) 

/9— >00 



The ground state solution that is positive and mono- 
tonically decaying is sought. One can show that 
< Qo{p;uj) < TT, therefore this solution is termed 
non-topological [7]. In contrast, there exist topolog- 
ical droplets that satisfy a related nonlinear eigen- 
value problem with the conditions 6o(0;a;) = tt, 
6o(p;w) 0. Their existence and stability subject 
to certain symmetric perturbations were studied 
rigorously in [21]. 

Assuming the existence of a localized solution, the 
following inequality holds 



< w < 1, 



(21) 



which can be seen by multiplying Eq. (20a) by p^ 

and integrating to find 



2co 



1 — cos0o) pdp - 



f 

Jo 



sin^ 60 pdp = 0. 



We immediately observe that 

_ /o°" sin^ Go p dp 



(22) 



u = 



2 /q°°(1 - cos Qo) pdp 



> 0. 



Using the fact that 2(1- cos 0) > sin^ 6 for < 
< TT, we also have oj < 1. 

The decay rate of Go can be obtained by expanding 
(20a) for large values of p and neglecting nonlinear 
terms which gives 

2 — ^j:r + (i - 1^) ©o ~ o . 



dp2 p dp 



givmg 



00 (p) = O 



Vp 



p > 1 . (23) 



From numerical computations, see, e.g., [10], it is 
known that the amplitude of the droplet increases 



with decreasing frequency. As the frequency ap- 
proaches zero, the droplet resembles a circular do- 
main wall with a transition layer that moves further 
away from the droplet center. The nonlinear, or- 
bital stability of nontopological, three-dimensional 
solitary waves in the context of anisotropic, two sub- 
lattice ferritcs were studied by use of a Lyapunov 
functional whose positive definiteness is guaranteed 
when [22] 

which is commonly referred to as a slope or VK con- 
dition after Vakhitov-Kolokolov's work on Nonlinear 
Schrodinger solitary waves [23]. Since N'{uj) is the 
Jacobian of the map from the physical variable to to 
the conserved quantity A/", we refer to (24) as a Ja- 
cobian condition as this will naturally generalize to 
the stability criterion for propagating droplets. The 
.Jacobian condition (24) has been numerically veri- 
fied for 2D stationary droplets [7], suggesting their 
orbital stability. Their robustness to perturbations 
in 2D, time-dependent numerical simulations [8] are 
further evidence of their stability. 

3. Propagating Droplet: Problem 
Formulation and Basic Results 

A two-parameter family of moving droplet solu- 
tions of (15) is considered in the form 

Q{x,t) = e{$,;uj,V), (25a) 
^{x,t)=u!t + (j>{ii;LO,V) , (25b) 

i = x-Vt, (25c) 



with 



where w is the droplet 's rest frequency and V is 
the droplet 's velocity. Insertion of this ansatz into 
Eq. (15) results in 

sin6» (w- V-V0) = 

- A6l-Fsin6' cos6» (l [V^l^) , (26a) 



sine V • V0 = V • (sin^eV^ 



(26b) 



The propagating nontopological droplet that we seek 
satisfies the boundary conditions 

lim 6»(C;w,F) = 0, lim V0($; w, V) = const , 

l€Koo |£|->-oo 

(27) 

with 9{^;(jj, V) < -K for all ^ e M^. Assuming decay 
for \^\ sufficiently large, 6* <C 1 and V0 ~ C/ a con- 
stant vector, Eq. (26b) is asymptotically 

-V ■ve = 2u -ve. 
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Therefore, 



lim Vm^,V) 



V 



(28) 



Due to the 5*0(2) invariance of the LL equation 
(10), it is sufficient to consider 



V = Vei, V>0. 



(29) 



For the rest of this work, expression (29) is assumed 
to hold. 

By use of (16), we observe that (26) is equivalent 
to the following 



60 

5_ 

J4> 



{ujM+VVi) 
{wJ\f+VVi) 



58^ 
58 



(30a) 
(30b) 



where V\ is the nonvanishing component of V when 

condition (29) holds. Equations (30) arc the neces- 
sary conditions for a minimizcr of the functional 



= £[e,(f)\-ujU[e,(j)\-vri[e,4>], (3i) 



where u> and V are Lagrange multipliers ensuring the 
constraints of fixed spin density M and momentum 
Vi - Propagating, localized structures in 2D isotropic 
[9] and 3D anisotropic [24] ferromagnets were studied 
numerically by minimizing a similar functional. Note 
that, 

While the variational formulation of the droplet leads 
to a parameterization in terms of A/" and V\ , for phys- 
ical applications, it is natural to parameterize the 
droplet in terms of the physical variables w and V. 
Solving Eqs. (26) for fixed to and V provides the map 
(w, V) i-> {Af, Vi). Invertibility of this map requires 
that the Jacobian 

dVi dN dAf dPi 



J 



(33) 



dV du) dV dto ' 

be nonzero. In the context of other magnetic and 
nonlinear wave models, the orbital stability of two- 
parameter solitary waves is determined by the Jaco- 
bian condition [25,26] 



J < 0. 



(34) 



Note that the Jacobian condition (34) degenerates to 
the Jacobian condition (24) for stationary droplets 
when V = because d^Pi = and dyVi > 0. 
Droplet stability will be discussed further in Sees. 4.1 
and 5.3. 

If 6^ (^) , 0* (^) are a minimizer of (31) , then dilation 
of ^ by the factor A gives 



A necessary condition for 0^{^), to be a mini- 

mum of (31) is therefore 



d'{l)= jsin^e^ Ai-VVi-2wJ\f = Q, 



(35) 



which can be derived more generally by direct manip- 
ulation of conservation laws for the Landau-Lifshitz 
equation [27] . Similar identities were studied by Der- 
rick [28] and Pohozacv [29] . Wc will make use of this 
relation in order to verify the validity of our numeri- 
cal computations in the Appendix. 

Another formulation of the boundary value prob- 
lem (26), (27), (28) can be derived by considering 
Eq. (18) with a traveling wave solution in the form 

w{xi,X2,t)= f{xi-Vt,X2)e''^^ 

= m,^2)e'^K (36) 

Inserting the ansatz (36) into Eq. (18) leads to the 
following nonlinear eigenvalue problem 



2r(/6 + 4) ,1-1/1^ 



1 + l/p ■'1 + l/p 

with / ^ as 1^1 — >■ DO. Letting 

/(Ci,6)=^^(6,6)e-^^«S 
then (37) becomes 

2U*[(«5, -Z|zi)2+«2J ^_ 



0, (37) 



1+ «• 



— u 



l + \u\- 



= 0. 



(38) 

Wc observe that for V ^ 0, u \s generally complex- 
valued. Assuming localization and linearizing 
Eq. (38) for |^| — >• oo, we obtain 



w ■ 



- 1 U = 0. 



Therefore, a necessary condition for localization is 



z/2(V,a;)>0, v{V,w) 



1-w- 



2 \ 1/2 



, (39) 



and the exponential decay rate of u is v{V,(jj). We 
discuss the consequences of this observation in the 
next section. 

3.1. The nonlinear dispersion relation 

From the boundary condition in (28), we observe 
that the propagating droplet is modulated by a spin 
wave background with wave number 

fc = J , (40) 



6 



wnl(A:, w) = w ■ 



(41) 



and frequency in the laboratory frame 

-=. + 2, 

In contrast to the hnear dispersion relation for ex- 
change spin waves in (14), the nonlinear droplet dis- 
persion relation in (41) depends on both the rest fre- 
quency w and the speed V . The condition (39) cor- 
responds to a negative nonlinear frequency shift 

WL(fc) — WNL(fc5 w) = 1 — OJ — fc^ = ^^(V, w) > 0. 

This yields an upper bound on the droplet's speed 
for a given rest frequency 

< F < V^^{ij) = 2Vl-a;, 

or, equivalent ly, an upper bound on the rest fre- 
quency for fixed speed 

4 



w < a;max(F) = 1 



y > 0. 



(42) 



In the limit w — 1, condition (42), via (40), re- 
quires that fc — > 0. The nonlinear dispersion relation 
(41) then coincides with the linear dispersion relation 
(14) for a spatially uniform, small amplitude excita- 
tion at the bottom of the spin wave band. Alterna- 
tively, in the limit k — > fcmax = Vl ~ we observe 
that ujNL{k,oj) 2 — uj = ojLikmax) corresponding 
to a propagating exchange spin wave. The speed V 
of the moving droplet is the group velocity Vg of the 
background spin wave, Vg = dtui^/dk = 2k = V. 
The phase speed of the background spin wave cjl /fc = 
uj/k + 2k can be greater than (when w > 0) or less 
than (when cj < 0) the droplet traveling wave speed 
V = 2k. In this work, we focus on the case uj > 0. 

The dispersion relations are illustrated in the lab- 
oratory frame and the comoving droplet frame in 
Fig. 2. The shaded regions correspond to the local- 
ization condition (39) hence localized, propagating 
droplets can be interpreted as nonlinear bound states 
lying below the linear spin wave band [7]. 

3.2. ID Droplet 

The nonlinear dispersion relation in Eq. (41) for 
a 2D propagating droplet is the same for a propa- 
gating ID droplet [7]. As will be shown, a number 
of 2D droplet properties are qualitatively similar to 
the ID droplet properties. Therefore, in this section 
we briefly review some basic results about moving 
droplet solutions of the LL equation in (1+1) dimen- 
sions. 

The Landau-Lifshitz Eq. (10) in one dimension, 
A — > dxx , is solvable by means of the inverse scatter- 
ing transform (see, e.g., [30] and references therein) 



« at 











t 





I 



■3.0 
P 



Wave number (k) 



,y=0 



Speed (V) 



Fig. 2. Dispersion relations in the laboratory frame (left) 
and the comoving droplet frame (right) with wavenumber 
k = y/2. The shaded regions represent a<;cessible droplet 
solutions. 

and admits soliton solutions. The single droplet soli- 
ton solution can be written as follows (see [12]) 



^ y2 _ 2 + 3w + VV^ + w2 cosh[2i/e] ^ 
cos © = , — — (43a) 



2-L0 + VW+^ cosh[2zy^] 
2 



+ tan 



4i^tanh(t/^/2) 



, (43b) 



where ^ = x — V t and v is given in (39). Thus, the 
ID soliton solution exists for all V > 0, cv satisfying 
(39). When V = 0, the rest frequency is limited to the 
unit interval (21). For u < 0,V = 0, the soliton (43) 
exhibits a phase singularity at the origin correspond- 
ing to topological behavior in this one dimensional 
setting. However, for V > 0, the necessary condition 
w < iOmax{y) in (42) has no lower bound. Therefore, 
nontopological solitons with arbitrarily negative rest 
frequencies exist even with vanishingly small prop- 
agation speeds. The phase gradient (43b) has the 
asymptotics 



$5(0) 



-V 



L0 + \cu\ + Vyi2\Lu\)'' 



< 



V 



«1, 



hence can be quite large when u; < 0. Our numerical 
method for computing 2D propagating droplets fails 
for modes with large momentum due to large phase 
gradients (see Sec. 5.2) so we are limited in exploring 
this regime of paramciter space in the two dimensional 
case. 

Explicit expressions for the energy £^^^\ the spin 
density A/'^^'^-' and the momentum 7"^^^' , as functions 
of V and w, for a single moving droplet solution, can 
be obtained from solution (43) via Eqs. (11), (12) and 
(13): 
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=2V4-4w-y2^ 
Af^^^^ = 2 tanh- 

p(iD)=^_2tan-i , , 

yV'i-4oj-v^, 

We observe the notable fact that the momentum 
is bounded, < V^^^'> < 27r, but the energy 
and spin density are not. 

4. Propagating Droplet: Asymptotic Results 

In this section, asymptotic methods are employed 
to study the behavior of the propagating droplet in 
the small amplitude and small speed regimes. 

4.1. Weakly nonlinear regime 

Small amplitude, propagating droplets were stud- 
ied asymptotically in [11]. In this section, we briefly 
rciview the results with a view toward their numerical 
verification in Sec. 5. 

Using a slowly varying envelope, small amplitude 
approximation the authors in [11] construct a propa- 
gating droplet bifurcating from the set of linear waves 
where the small parameter is the deviation from the 
linear band edge (39) 

< iy(V,uj)'^ < 1. 

The droplet solution for Eq. (18) takes the approxi- 
mate form 



w{xi,X2,t) 



-i[V{xi-Vt)/2-u>t] 



V2(2 - w) (44^ 

with > having the meaning of the inverse half 
width. As the band edge is approached, i/ — >■ 0, 
the propagating droplet is approximately a wide, 
small amplitude, propagating Towncs mode [31], the 
unique, positive, monotonically decaying solution of 

—u'-^——Urp — + Ut = , 



R 

u'j,{0) = 0, 



lim Ut{R) = 0. 



Using (44) and its next order correction, the func- 
tional (31) can be expanded as [11] 



where 



1 + V^/A 



2 ^ 4 



Wt\\\2( 



Urp ttrp 



2(1 + 1/2/4)2 

11.7009, 
■RdR 



(45) 



1.35825. 



Combining (45) with the relations (32), the propa- 
gating droplet's spin density, momentum, and energy 
in the weakly nonlinear regime are approximately 



N^ 



Nt 



1 + ^2/4 
NtV 
2(1 + 1/2/4) 



1 + 



,6(i-y2/4) 



(1 



T/2/4)2_ 
l-y2/4 

'(1 + 1/2/4)2 

1 



S^Nt 



J2b{l 



i + y2/4 

1/2/4) 



V (1 - 1^2/4)3 

1 + 3^2/4 



(46a) 
(46b) 

(46c) 



+ 



1 



(1 + ^2/4)2 i + V-2/4^ 

These results are numerically verified in Sec. 5.2. 

For the initial value problem in the weakly non- 
linear regime, one can derive at leading order the 
time-dependent cubic Nonlinear Schrodinger equa- 
tion (NLS) from Eq. (18) for the slowly varying en- 
velope of a carrier wave. In two spatial dimensions, 
the NLS equation exhibits critical blow up with the 
Townes mode playing an important role [32] . For ini- 
tial data with L?{E?) norm less than |]ut||2, the so- 
lution disperses to zero as i — > oo. In other cases, 
e.g. for data close to ut in L'^{M?) but with larger 
norm, the solution blows up in finite time. Thus, the 
Townes mode is unstable to small perturbations. It is 
natural then to inquire into the stability of the prop- 
agating droplet for the LL Eq. (18). The Jacobian 
condition (34) was verified in [11]. Insertion of (46a) 
and (46b) into Eq. (33), leads to 

""^^ •Oiu^), (47) 



J = 



2(2 



suggesting that the higher order corrections to 
the NLS approximation stabilize the propagating 
droplet. Further discussion of stability in the large 
amplitude case is presented in Sec. 5.3. 

4.2. Small propagation speed 

We now consider Eqs. (26a) and (26b) with bound- 
ary conditions (27) and (28) for small values of the 
droplet speed Q <V ^1. Polar coordinates are used 
(see Fig. 1) 

S^i ^ xi — Vt — p cos , £,2 = X2 = p s^n^p . (48) 
A solution is sought in the form 
^ = eo(p;w) + 0(l/2), 

(/) = V (^$i(p,^)-^pcos^^ +0(y2)^ i"^^) 

where 8o is the solution of Eq. (20) for the stationary 
droplet with < w < 1 and and (p are given in 
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the ansatz (25). The aziniuthal angle (p in Eq. (49) 
should not be confused with the angle (f of Eq. (48); 
see Fig. 1. 

At first order in V, the function $i satisfies 

sinGo A$i+2 9g cosSq = (cos6o-l)6o cost^, 

(50) 

subject to the boundary condition 
lim V$i = 0. 

p—hoo 

Here and hereafter the prime indicates derivation 
with respect to p. Equation (50) admits separation 
of variables in the form 



*i(A</5) 



$l(/9) COSt^ 



where $i is the solution of the ordinary differential 
equation (ODE) 

l>i = (1 - cos 6o)6o . 
d^ 



L --—--— + — 
* dp2 p dp /)2 



V^ip), 



V$(p) = — 6q + (cos 6o — w) cos 00 , 
subject to 
d 



(51a) 

(51b) 
(51c) 



lim 

P^oo dp 



y^ePVi— $^(p) =0, $i(0)=0, 



(51d) 

from the asyniptotics of the stationary droplet (23). 

Equation (51a) is solvable for any choice of the 
parameter < a; < 1. This can be seen as follows. 
Standard arguments {e.g. using variational methods) 
show that L$ is a positive definite operator. The so- 
lution to (51a), 

$1 = L^i{(i-coseo)e;,} , 

satisfies the boundary condition (51d) as shown by 
the asymptotics, 

$1 = O (p) < p «; 1 . 

Therefore, from the ansatz (25) and (49), we find that 
an approximate droplet solution for small V consists 
of the stationary solution modulated by the spatially 
varying phase 

5. Propagating Droplet: Numerical Results 

In this section, a numerical method is introduced 
to solve Eq. (37) for propagating droplets. The so- 
lutions' properties are investigated and its stability 



verified by the analytical criterion J > and time 
dependent numerical simulations. 

5.1. Iterative Method to Compute Propagating 
Droplets 

A solution to the nonlinear eigenvalue problem (37) 

is computed by the spectral rcnormalization method 
[33] , a generalization of the Petviashvili method [34] . 

We begin the procedure by rephrasing Eq. (37) as 
follows 

ujf + iVU,+Af = N[f], (52a) 

with 

2r(/x^+/x'j+/(i-i/n 



m 



(52b) 



Fn+l — 



1 + 

and, upon taking its Fourier transform, results in 

(co-Vki-k^) f = N[f], 

where / = J^[/](k) is the Fourier transform of / and 
K = {ki,K2) is the Fourier wavenumber, with |k.| = 
K. The iteration 

{uj-Vk,-k^-c) = N[f^] - cfn , (53) 

is introduced so that 

LU — C — V Kl — ' 

where c is a constant real parameter, chosen in order 
to avoid divergence of the iteration - its value is ex- 
perimentally obtained, see [33]. The value of c affects 
the convergence speed of the iteration and must be 
taken larger as the droplet amplitude is increased. In 
most of our computations, we find c = 10 to suffice. 

In order to iterate the sequence {/n}, we scale Fn+i 
according to 

fn+l = Vn+l Fn+l , (54) 

and choose the complex parameter r]n+i so that the 
projection of Eq. (52) onto fn+i is satisfied. Insert- 
ing (54) into Eq. (52) and taking the I/^(M^) inner 
product with fn+i results in 

dFn+i 



I 



Fn 



iVF 



n+l 



dx 



(55) 



|VF„. 



dxi dx2 



Vn+l\ X 
2" 



(%^) +(S%r) 



+ 



I 



l + |r/„+i|2|i^„+i|2 

l+|7?„+lP|F, 



dxi dx2 



\Fn+l\ 
12 



da;i dx2 



n+l\ 



Equation (55) determines [r^n+i). 
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The numerical computations are implemented on 
the finite, discrete grid 



a;i,2,- = -i^ + (j - 1) A , j = l,2, 



with 



A = 



2L 



Due to the localization of -F„+i, derivatives are esti- 
mated using the fast Fourier transform. We take L = 
50 and A'^ = 1536, implying A « 0.065, imless other- 
wise stated. To obtain droplets with large speed, we 
use continuation in V. For small V , the initial guess 
/o is the stationary droplet discussed in Sec. 2.2 with 
the same ui. For large V. /o is a previously computed 
propagating droplet with the same oj but smaller V. 
The integrals in (55) are computed using the spec- 
trally accurate trapezoidal rule. 

Equation (55) is solved for |/7„+i|^ using Newton's 
method and the phase of r]n+i is determined by the 
condition that /„+i have a vanishing phase at the 
origin so that 

/n+i = \Vn+i I K+1 e-' [^--s , (56) 

This implies that the projection of the magnetization 
at the origin onto the plane points along ei . 
It is helpful to introduce the residual 

r[fn] = [^n - - - K^) f 

which is zero at a droplet solution. We iterate from 
/o until 



Wfn- 

\\r[fn 

\Vn - 



/n-l||oo < tol/ - 
II oo tol J- , 
77„_l| < tol^ , 



where we set tol/ = tol^ = toV = 10~^, unless oth- 
erwise stated. 

Further details and validation of this numerical 
scheme are presented in the Appendix. 

5.2. Droplet Properties 

Using the numerical scheme presented in the previ- 
ous section, we computed a large number of solutions 
to the boundary value problem (37) for a propagat- 
ing droplet with parameters lying in the set 

S= {{V,uj) I < y < 2\/r^ - 0.25, 
0.275 <u< 0.95}. 

This representative set of droplet parameters was 
chosen because accurate solutions were obtained 
without having to change the numerical parameters 
of our scheme. Outside of this set, either higher reso- 
lution or larger domains were often required to accu- 
rately resolve the droplets (see Appendix). Droplets 



with w < were computed and exhibit similar prop- 
erties to those in S. We were unable to compute 
large momentum droplets for particular V with oj 
sufficiently small or negative due to inherently large 
phase gradients. Illustrative example droplet solu- 
tions are shown in Fig. 3. As oj and V are increased, 
the droplet amplitude, defined as 1 — min (m^) ~ 
1 — ma (0,0), decreases. The spin wave background 
is evident, exhibiting a decreasing wavelength as V 
is increased or as lu is decreased. This is in accord 
with the nonlinear dispersion relation discussed in 
Sec. 3.1. The precessional motion of a propagating 
droplet is depicted in Fig. 4 in the comoving frame. 

Using Eq. (19), the energy £, momentum Pi, and 
spin density Af were computed over S. Contour plots 
of these quantities arc shown in Figs. 5(a)-5(c). Ad- 
ditional values of the conserved quantities outside 
S were obtained by interpolation using the weakly 
nonlinear asymptotics (46) . Both the energy and the 
spin density are monotone decreasing functions of V 
and OJ. On the contrary. Fig. 5(b) shows that the mo- 
mentum dependence is not monotonic. For fixed val- 
ues of to, the momentum Vi exhibits a maximum for 
V > 0. This "ridgeline" is shown as the dashed curve 
in Fig. 5(b). The droplets with the largest momen- 
tum exhibit the largest phase gradient near the ori- 
gin. See, for example, the case V = 0.212, to = 0.3 in 
Figs. 3 and 4 which is on the ridgeline of large mo- 
mentum droplets. 

Figure 6 depicts the numerical verification of the 
weakly nonlinear results (46) for fixed w = 0.4. The 
results agree near the band edge where ^max — V <^ 
1 and we have verified that the error exhibits the 
expected i^'* scaling. 

The droplet's amplitude dependence is shown 
in Fig. 5(d). As observed earlier, larger amplitude 
droplets correspond to smaller speeds and rest fre- 
quencies. 

The droplets shown in Figs. 3 and 4 show negligible 

signs of asymmetry in their localization extent. To 
further investigate this, we define the droplet widths 
in the Xi and X2 directions according to the standard 
deviations 



widthi 



width2 



1 



x\ (1 — cos 6) dx 



x\ (1 — cos 9) dx 



(57a) 



(57b) 



A contour plot of the quantity widthi for (V, oj) G S 
is shown in Fig. 5(e), exhibiting non monotonic be- 
havior for increasing V and w, especially for large 
amplitude modes. The droplet's width diverges ac- 
cording to 0{l/i'{V,uj)) as the hnear spin wave band 
is approached < i/ <C 1 (see Eq. (44)). 
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Fig. 3. Computed droplet solutions for several choices of oj and V. The color contour plot represents the magnitude of ms and 
the arrows represent the projection of the magnetization on the plane. To facilitate visualization, the length of the projected 
field is normalized to the largest value in each frame. 
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Fig. 4. Precessional motion of a propagating droplet with = 0.3 and V = 0.212 in the comoving frame. Each frame depicts a 
different instant in time over one period 2tt/lo. 



The asymmetry, defined as (width2— widthi)/"widthi, 
is a nonnegative quantity as shown in Fig. 5(f). 
Even though the asymmetry does increase - ap- 
proximately with the momentum - it is very small 
for droplets with parameters in S, reaching a maxi- 
mum of approximately 0.066. Wide, small amplitude 
droplets are symmetric as shown in Eq. (44). 



We now comment briefly on the localized exci- 
tations of Equation (10) observed in [10]. These 
states were formed by starting with a stationary 
droplet initial condition /o(/o;cli) — sin Op (p; a;)/ [1 -I- 
cos0o(/9;w)] to Eq. (18) superimposed with a spa- 
tially varying phase in the form 
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(e) Width (xi-direction): widthi (f) Asymmetry ratio: 
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Fig. 5. Contour plots of the droplet's energy, momentum, spin density, amplitude, width in the x-direction and asymmetry 
ratio, computed from numerical solutions of the boundary value problem. 
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(a) Energy, £ 



(b) Momentum, Vi 
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Fig. 6. Comparison of the weakly nonlinear asymptotics (46) with numerical computation of £, V, and M for uj = 0.4, 

Vmax ~ 1.55. 



w{x,t^Q)^e'-'^-^h{\x\-u). (58) 

This initial condition to Eq. (18) was allowed to 
evolve in time and a moving, localized structure was 
observed. The frequency of precession w, speed V , 
and energy E' associated with several initial condi- 
tions reported in [10] are given in Table 1. For com- 
parison, the energies £ of directly computed propa- 
gating droplets with the same frequency and speed 
are also listed in Table 1. The droplet approximation 
(58) leads to overestimated energies which are less 
accurate for larger speeds. This behavior coincides 



with the small V asymptotics of Sec. 4.2 which show 
that the ansatz in Eq. (58) is approximately valid. 

5.3. Propagating Droplet Stability 

The nonlinear, orbital stability of droplets has been 
supported by verification of the Jacobian condition 
(24) in the stationary [7] and weakly nonlinear [11] 
(recall Eq. (47)) cases. In this section, we confirm the 
stability of large amplitude, propagating droplets by 
numerically demonstrating that J < and by per- 
forming time-dependent numerical simulations of the 
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U) 


S' 


S 


0.067 


0.50 


25.6 


25.5 


0.13 


0.49 


25.9 


25.8 


0.24 


0.47 


27.0 


26.9 


0.3 


0.46 


28.9 


27.8 


0.35 


0.44 


31.4 


28.8 



Table 1 

Comparison between the energies E' computed in [10] by 
solving the initial value problem (58) with the energies S of 
droplets computed in this work. 




Fig. 7. Numerical computation of the negative Jacobian — J 
in Eq. (33) demonstrating that the map {uj,V) i— (A/", "Pi) 
is invertible (J ^ 0) and the Jacobian condition J < for 
stability is satisfied. 

LL Equation (10) with initially perturbed propagat- 
ing droplets. 

Figure 7 exhibits the numerical evaluation of the 
negative Jacobian — J (33) for {V,uj) e S. Since J < 
0, we have numerically confirmed the invertibility of 
the map {lo,V) ^ {N ,Vi) and the Jacobian condi- 
tion (34) for droplet stability. The invertibility of the 
map is not obvious given the lack of monotonicity of 
Vi in Fig. 5(b). Note also that the Jacobian is mono- 
tonic in both uj and V within the set S and in the 
weakly nonlinear regime (Eq. (47)). 

Recalling that the Jacobian condition J < for or- 
bital stability of solitary waves has only been demon- 
strated in other nonlinear wave systems [25,22,26], 
we seek further evidence of stability by performing 
time-dependent numerical simulations of perturbed 
droplets. To solve Eq. (10) numerically, we use the 
method of lines with a pseudospectral, Fourier dis- 
cretization in space, the same as that used in our it- 
erative numerical scheme of Sec. 5.1. For time step- 
ping, we used an adaptive time stepping, second or- 
der Runge-Kutta ODE solver. Numerical parameters 



we used were L = 50, = 768, giving A « 0.13. Be- 
cause the time stepping method does not preserve the 
length of the magnetization vector, after each time 
step the solution is rescaled so that \m\ = 1 is pre- 
served. The initial conditions consist of a computed 
droplet perturbed by noise (explained below). In all 
computations, we have verified that the quantities 
(11), (12), and (13) change at most by 0.74% over 
the course of the simulation. 

The initial conditions are constructed as follows. 
Uncorrelated, mean zero, normal random variables 
with variance 0.2 are sampled at each grid point to 
construct a complex valued noise function n. A low 
pass filter is applied with cutoff wavenumber 6 yield- 
ing fi. Adding a droplet, /, computed as in Sec. 5.1 
to the noise yields the perturbed droplet f = f + n 
in stereographic form. Inverting the relation (17) 



1 



mi -I- im2 = (1 + "Is)/ , 



1 + 

yields the initial condition for the LL equation (10) 
in vector form. 

An example solution at i = 50 is shown in Fig. 
8. The perturbed droplet 's localization structure 
has been maintained. Furthermore, it's propagation 
speed and precessional coherence have also been 
approximately maintained. To demonstrate this, we 
show in Figs. 9(a) and 9(b) the position and phase, 
respectively, of the minimum of the computed so- 
lution as a function of time. For comparison, we 
have included the trajectories of these quantities for 
unperturbed droplets. Stability for the considered 
droplets is clearly achieved with only small changes 
in the propagation speed and rest frequency. These 




Fig. 8. Perturbed droplet with (V,uj) = (0.7,0.3) at t = 50. 
Compare with unperturbed case in Fig. 3. 

results suggest that droplets are robust, approxi- 
mately maintaining their precessional frequency and 
propagation speed, even after experiencing signifi- 
cant perturbations. 
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The droplet iterative scheme of Sec. 5.1 converges 
for parameter values satisfying i/^(V,a;) < 0. How- 
ever, as expected, such solutions are unstable due to 
linear spin wave resonances. 




10 20 30 40 50 

t 

(b) Phase at solution minimum. 

Fig. 9. Time evolution of the xi position (a) and the phase 
(b) at the minimum of three perturbed droplets (solid curves) 
compared with unperturbed droplet evolution (dashed lines) . 



6. Discussion 

In many respects, the 2D propagating nontopo- 
logical droplet considered here has much in common 
with the ID propagating droplet soliton reviewed in 
Sec. 3.2. Moreover, by rewriting Eqs. (26) in polar 
coordinates, it is immediate to observe that, for large 
values of the polar radius p, the droplet behaves as a 
1-dimensional object (namely, depending only on p, 
with the velocity V replaced by V costp). 

One interesting feature of the ID case is the devel- 
opment of topological structure, a phase singularity 



for V = and w < 0. While it is not possible for a 
single, topological solitary wave to stably propagate 
[18], the copropagation of pairs with opposite charge 
are possible [27]. The bifurcation of nontopological 
droplet solitary waves to pairs of topological struc- 
tures have been predicted in isotropic ferromagnets 
[9] and in 2D, uniaxial, easy-plane ferromagnets [27]. 
It is natural to ask whether there exist propagat- 
ing solutions with some local topological structure 
for the LL equation (10). Time-dependent numerical 
simulations of an appropriately chosen initial value 
problem have apparently yielded such structures [10] . 
Our computations with the stereographic Eq. (37) 
are limited to nonsingular, hence locally nontopo- 
logical solutions. Furthermore, the large momentum 
states we have computed exhibit large phase gradi- 
ents near the droplet center, suggesting the possibil- 
ity of a bifurcation into local topological structure for 
sufficiently small or negative rest frequencies. 

Another, necessarily multidimensional property of 
droplets was studied in [10] where ninety degree scat- 
tering was numerically observed for the head-on col- 
lision of two "accelerated" stationary droplets with 
the same frequency and opposite velocities. Merg- 
ing accompanied by spin wave radiation was also ob- 
served for certain initial configurations. The prop- 
agating droplets constructed in this work could be 
used as initial conditions for a direct study of 2D 
droplet interactions. 

7. Conclusions 

Through numerical and asymptotic means, we 
have demonstrated the existence and studied the 
properties of moving droplets across a wide range of 
precessional frequencies and speeds as well as their 
stability to perturbations. Along with the rise of 
nanomagnetism research and applications [5], this 
work encourages further study of localized magnetic 
structures including droplet propagation in the pres- 
ence of damping, external magnetic fields, nonlocal 
dipolar fields, and injected spin current. 

Appendix: Validation of the Numerical 
Scheme 

We validate the numerical scheme described in Sec. 
5.1 in several ways. First, the Derrick identity (35) 
was checked for the computed droplet modes with 
(y, oj) e S. The largest deviation from zero of the 
quantity dl (1) in (35) was found to be approximately 
5.3 • 10^^, which occurred for the stationary droplet 
y = 0, w = 0.9, likely due to the finite domain size. 



14 




Fig. 10. Absolute value of the Fourier coefficients, normalized 
by the maximum, for a single droplet with = 0.3 and 

V = 0.212. 

Second, we verify that the resolution of the em- 
ployed discretization is sufficiently fine to resolve the 
solution's features. As can be observed in Fig. 10, we 
confirm that the relative magnitudes of the Fourier 
coefficients decay to machine precision for a large mo- 
mentum droplet solution with V ~ 0.212 and uj ~ 
0.3. We find that all droplets with {V,uj) E S show 
similar behavior with a large momentum droplet ex- 
hibiting the slowest decay to at most 4.9 • 10~^ for 
large k, still an acceptable value. This suggests neg- 
ligible aliasing effects so that we have properly re- 
solved the solution. 

Finally, to determine a domain size and grid reso- 
lution that accurately resolves the droplets, we com- 
puted a very accurate, yardstick solution, /yg, with 

V — 0.7 and uj — 0.3. The numerical parameters 
used were L = 50, iV = 2048, c = 10, and tol/ = 
to\r = tol,, = 10"^^ so that A « 0.0488. We compare 
fys with the same droplet mode obtained on different 
grids (A,L). The relative error - /||oo/||/ysl!oo 
is computed using sine interpolation [35] of the yard- 
stick to the grid for / and plotted in Fig. 11. Domain 
truncation error decreases with increasing domain 
size until discretization error dominates. The small- 
est grid spacing A = 0.075 with i > 40 yields rel- 
ative errors below the precision of the yardstick tol- 
erances. The iterative computations of droplets with 
{V,uj) G S employ the conservative values (A,L) = 
(0.065, 50), accounting for the larger phase gradients 
that occur. 
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